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Abstract. We present a family of closed form inversion formulas in ther- 
moacoustic tomography in the case of a constant sound speed. The formulas 
are presented in both time-domain and frequency-domain versions. As special 
cases, they imply most of the previously known filtered backprojection type 
formulas. 

1. Introduction and statement of main results 

Thermoacoustic tomography (TAT) has recently attracted considerable attention 
as a promising method of biomedical imaging [28l EU [201 ED US] • We give here a 
brief introduction to the mathematical model of TAT. 

An electromagnetic (EM) pulse in visible light or radiofrequency range is sent 
through a biological object of interest. A fraction of the EM energy is absorbed by 
the tissues throughout the object. It is known that cancerous tissues absorb much 
more EM energy than healthy ones (e.g., [2H1II6]). Thus, knowing the absorption 
distribution function a{x) would be useful for detecting cancerous tissues. The EM 
energy absorbed by the tissue causes its thermoelastic expansion, which in turn 
releases a pressure (ultrasound) wave u{x,t) propagating through the body. This 
pressure is then measured by transducers located on an observation surface S 
surrounding the object. The initial pressure f{x) = M(a;,0) is roughly proportional 
to a{x). One now concentrates on the recovery of f{x) from the measured data 
g :— u\sx[o,oo)- There has been a number of papers devoted to this reconstruction 
(see reviews in [131 HH 113 E])- AH known closed form inversion formulas assume 
that the object is acoustically homogeneous, i.e. the sound speed is constant. We 
will assume that, after choosing appropriate units, the sound speed is equal to one. 

The standard mathematical model of TAT is [26l [71 [22] 

{uttix, t) - Au{x, t) = 0, xe IR", t > 0, 
uix,0)^f{x), ut{x,0)^0, 
u{y,t)= g{y,t), for y G 5, t>Q. 

Here g{x, t) is the measured data and / is the function to be reconstructed. In 
other words, one needs to invert the operator T : f i—f g. Since T is known to 
be invertible from the left only, different inversion formulas exist. In this text, we 
develop a family of explicit closed form inversion formulas in the case when S is 
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the unit sphere centered at the origin and / G C^{B), where B is the open unit 
baU enclosed by S. 

We will also need to deal with some other operators closely related to T. Consider 
the wave equation 

{uu{x, t) - Ait(x, t) = 0, a; G H", t > 0, 
w(x,0) = 0, ut{x,0) ^ fix), 
u{y,t) = g{y,t), for y G 5, t > 0. 

One can now define the the operator V that maps function / into g: V{f) = g. 
Then V is related to T as follows: 

t 

V{f){y,t) = j Tif){y,T)dT and r(/)(t) = dtVim). 



We also introduce the spherical Radon transform TZs with centers on S by 
7^s(/)(y,^)= J f{y + ruj)r"-^d<j{uj), ioYS.\l{y,t)eSxlR+. 

Finally, Als is the spherical means operator with centers on S: 

Ms{f)iy,t)^— f /(y + r^)daH, foraU(y,t)e5x]R+. 

In these formulas, da{uj) is the standard surface measure on the unit sphere S"^^ C 
IR" and ujn is the total measure of the sphere. 

In various works cited below, different operators from this list were considered. 
However, due to known exphcit connections between T and V, TZs, and Ms (e.g., 
[SI US]) J inversion formulas for these operators are closely related. 

The first such formulas, in odd dimensions, were obtained by Finch, Patch and 
Rakesh [TU] using some trace identities for the wave operator. Xu and Wang [57] 
derived a different formula for n = 3 by working in the frequency domain. A formula 
for all dimensions, which coincides with that of [27j when n = 3, was presented by 
Kunyansky [18]. Its derivation is based upon some symmetry relation for special 
functions. By applying the same method as in |T0], Finch, Haltmeier and Rakesh [9] 
obtained inversion formulas for even dimensions, which involves the data measured 
for an infinite time period. The authors of 9J also derived another type of inversion 
formulas for even n, which uses only the data measured for a finite period of time 
(which we will refer as " finite-time formulas for even dimensions" ) . 

In spite of availability of several types of closed form inversion formulas, some 
questions have remained unanswered. For instance, it was not clear, what is the 
relation, if any, between formulas of [HI [TO] and [TOl dZ], (which are known to be 
not equivalent outside the range of the operator T). The same applies to the two 
types of inversion formulas derived in '9] for even dimensions. 

The goal of this work is to obtain a unified family of closed form inversion 
formulas, which in particular would produce formulas of [lOl [9l [HI [21] as particular 
cases. 

We now formulate the main results of the paper. Consider the function 

Tl-2 

(3) «(^'^) = i(2^) ^ 
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where iji^^a is the Hankel function of the first kind. Then (e.g., pj) ^{x,y,X) = 

2 

G{\x — y\, X) is the solution for the Helmhohz equation 
(4) AUiy) + X^U{y)^-6iy-x), 

obtained by limiting absorption. For odd n and any s > 0, G{s, .) G C°°(IR). For 
even n, G(s, .) e C°°(IR\ {0}) with a logarithmic singularity at zero [1]. 

Let g{y,t) be the function in ([T|) that represents the thermoacoustic data and 
go{y, t) be its even extension with respect to t. Let also go be the Fourier transform 
of go with respect to time: 

30(2/, A) = ^J 9o{y,t)e'^'dt. 
u 

Then, our family of inversion formula reads as follows: 

Theorem 1.1. Suppose that f S C^{B) and g = Then for any x £ B and 

^ e K", the following equality holds 



(5) fix) = -2 / I / Gis,X)go{y,X)dX 



ds 

ii 

where G is the complex conjugate ofG. 



^y-'^'-^^daiy), 



\x - y\ 

= \x-y\ 



One can rewrite this inversion formula without going to the frequency domain. 
Namely, we define a transform W of a function v G C°°[0,oo) as follows (as long 
as the expressions involved make sense): 

cni '-fj ^ I-^dt, if n is even. 




J , if n is odd. 



(-i)~ 



H— , if n is even, 



„^_i , if n is odd. 

2(27r)T~ 



(As noted in Section [71 W is a known transform that intertwines the second deriv- 
ative and the Bessel operator.) One can now obtain another representation of the 
inner integral in ([5]): 

''G{s,X)goiy,X)dX = Wig){y,s), 



H 

where W(g)(y, s) :— W{gy){s) with gy{t) := g{y, t). Then Theorcm ll.ll is equivalent 
to 

Theorem 1.2. Suppose that f G C^{B) and g — T{f). Then for any x £ B and 

^ € H", the following equality holds 



(8) fix) = -2 / {±W{g)iy,s) 



{y~x,y~^) 

1 1 — d<^{y)- 
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One can look at the inner integral in ([5|) in a different manner yet. Namely, due 
to the relation between T and TZs^ it can be shown that 

oo oo 

j G(s,A)go(y,A)dA= ^ j XTZ{s,X) J TZs{f){y,r)I{r, X)dr, 

R 

where TZ and 2 are the real and imaginary parts of G. Then Theorem I f . 1 1 implies 

Theorem 1.3. Suppose that f e C^{B) and g = TZs{f)- Then for any x £ B and 

^ € H", the following equality holds 



s 



whe 



Kn{y,s) = / A7^(s,A) / giy,r)I{r,X)drdX 



The key ingredient in the proof of Theorem 11.11 is the following identity, which 
is equivalent to a range description of operator T (see Remark [ 



Theorem 1.4. Suppose that Q is a ball, f e C^(f7), and u solves {Ip. Let uq be 

the even extension of u with respect to t. Then for all x ^ ^ we have 



(10) j J G{\x-ylX)uoiy,X)dXdaiy) =0. 

on R 

This statement implies that adding to any inversion formula an expression 



Aj j G{\x-y\,X)goiy,X)dXda{y) 

S R 

with an arbitrary operator A also produces an inversion formula (since the added 
term vanishes on the range of the operator to be inverted). In particular, if A 
is the operator of multiplication by an arbitrary function (/3(x), one can add the 
expression 



^(x) J J G{\x - y\,X)goiy,X)dXda{y). 
S R 

to (HI and ^ to obtain the following inversion formulas: 

Corrolary 1.5. Suppose that f e C^{B) and g — T{f). Let (p be an arbitrary 
function defined on B. Then for any x d B and ^ G IR", the following equalities 
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hold 



(11) + ^{x) / K,,{y,\x^y\)da{y) 



s 



fix) = -2 / l^wigMs) 



s=\x-y\ \x y\ 



(12) + ^{x) J W(g)iy,\x-y\)daiy). 

s 

Now, by suitable choices of ^ and (p in pT|) and ([H]), one can recover the inversion 
formulas known in the literature. If we let S, = x, and if = — 2(n — 2), —2{n — 1), 
or —2n in (|12p . we obtain correspondingly the formulas derived in [HI Theorem 4] 
and [lOl Theorem 3]: 

Proposition 1.6. Let f G C^{B). Then for any x G B: 

(13) fix) = -2iV*td^Vf)ix) 

(14) fix) = -2{V*dttdtVf)ix) 

(15) fix) = -2iV*dftVf)ix) 

Here dt = is the derivative with respect to t, and V* is the L"^— adjoint of V . 

Choosing ^ = x and if = Q in (jlip . one gets a finite-time inversion formula for 
even dimensions similar to the second one in [9| Theorem 2]: 

Proposition 1.7. Assume that n is even, f G C^iB), and g = Msif)- Then for 
all X G B , fix) is equal to 



(-1)"2 UJrt 

(27r)" 



2 r 1 

-j^ \ n — 1 



drr dr- r'^-'g 



iy, r) In — \x — drda{y), 



S 



where LUn is the surface area of the unit sphere S^' 

In particular, when n = 2, we obtain the second formula in [51 Theorem 1]: 

Corrolary 1.8. Assume that n = 2, / e C^{B), and g — Msif)- Then for any 
X e B, 

2 

(16) = J[drrdrg]iy,r)lny-\x-y\^\drda{y). 

S 

Finally, let 

(S) Nr^ (S) 

Jis) = ^E^, Nis) = ^^, 

S 2 S 2 



n— 2 

:uiiiaiiii iuiic lions ui oiuei 

ing ^ = and = in pT|). one arrives at: 



where J n-2 and A^.1-2 are the Bessel and Neumann functions of order 2—1. Choos- 

00 ^ 



6 



LINH V. NGUYEN 



Proposition 1.9. Assume that f G C^{B) and g — TZs{f)- Then for any x ^ B, 

-1 



(17) f{x) = r,,o^n~i ^x I n{y)kn{y,\x-y\)da{y). 



S 

where 

oo 2 

fc„(j/,s) = j A2"-37V(sA) j g{y,r)J{r\)drd\. 



This inversion formula is equivalent to the one obtained in [18'. Indeed, let 

oo 2 

h{y,s) = j \^^-^N{s\) j g{y,r)J{rX)drd\ 



oo 2 

A2"-3j(sA) j g{y,r)N{r\)drd\. 



It can be shown that h[y,s) — 2k„{y,s), and thus Proposition 11.91 implies the 
following result of [TSj : 

Proposition 1.10. [TB] Assume that f G C^{B) and g = TZsif)- Then for any 
X & B, 

(18) fix) = 4(2~)Li V. I n{y)h{y, \x - y\)da{y). 

s 

Remark 1.11. Formulas of fTOllllllHlET] do not reconstruct f{x) inside S correctly 
if a part of support of f lies outside S. This feature is absent in other reconstruction 
methods such as time reversal (the readers are referred to |171 for these discus- 
sion). I would be interesting to see whether formulas contained in Corollary \1.5\ are 
any different in this regard. 

The paper is organized as follows. In Section [2l we derive Theorem 11.11 from 
Theorem 11.41 In Section [31 we show that Theorem 11.21 equivalent to Theorem 11.11 
and Theorem 11.31 is implied by Theorem 11.11 Two proofs of Theorem 11.41 are 
presented in Section [D In Section O Propositions II. 6[ II. 7[ II. 91 and II. lOl are derived 
from Corollary 11.51 using the aforementioned choices of ^ and (p. Proofs of some 
auxiliary results are given in Section [6l Finally, some remarks are provided in 
Section [71 

2. Derivation of Theorem 11.11 from Theorem 11.41 

Let u solves ([1} and uq be its even extension with respect to t. Then, uq solves 
the wave equation on IR" x IR,. Therefore, 

(19) A2Mo(2/,A) + AMo(y,A) = 0. 
Due to ([ll and (fT9|) . we have the Green's identity: 



(20|to(2:,A) = - 



dG{\x-y\,X) duQ{y,X)-^ ' 
7^ uo(y, A) G{\x - y|. A) 



S 



d(j{y). 
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Since f{x) = uo{x,0), the Fourier inversion gives 

f{x) = / uq{x, X)dX. 



Due to (j20|), we get 
(21) 



dGi\x-ylX) 
dv,, 



uoiy, X)dXda{y) 



S H 



^^fl^Gi\x-ylX)dXdaiy). 



S IR 



Let S'(0, r) be the sphere centered at the origin with radius r. For any function 
H e C^(IR"), by changing variables, we derive 



d 
dr 



I 



ii{y)do{y) 



V 5(0,r) 



d 

dr 



- / H{ry)da{y) 



That is, 



H{y)da{y) 



Applying this equality for 



dH 
dv,, 



{y)da{y). 



H{y) = / G(|y-x|,A)uo(y,A)dA, 



we get 



du,, 



G{\y - x\,X)uo{y, X)dXd(j{y) = 



IR, 



d 
dr 



Gi\y - x\,X)uoiy, X)dXda{y) 



S{0,r) u 

Due to Theorem II .41 the right hand side is zero. Thus, 

^ / G{\y - x\,X)uo{y, A)dAda(y) = 0. 



du, 



IR 



Therefore, 



(22) 



dGi\y-xlX) 



S IR 



uo(y, X)dXda{y) 



s u 



Gi\y-xlX)^iy,X)dXdaiy). 
avy 
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Combining this and ([2T|) . we get 



ff \ oil dG{\x-y\,\) , 
/(x) = -211 — uo{y,X) 



S TR. 



(23) = _2 / / G,(|x-y|,A)Ko(y,A)^^— 

\x-y\ 



s u 



Here Gs(s, A) is the derivative of G(s, A) with respect to s. 
Applying Theorem II. 41 for fl = B, we get 



(24) J J G(|x-y|,A)iio(2/,A) = 0. 

s u 

Taking the derivative with respect to x of the above identity along the direction ^, 
we obtain 

G.(|x-y|,A)iio(2/,A)fc^=0. 

\x - y\ 

S B, 

Subtracting this equality from (|23p . we conclude that 

fix) = -2 / /G.(|a;-y|,A)«o(y,A)fc^i^-^dAda(y) 

\x-y\ 



S H 



= -'^ji^J G{s,\)go{y,X)dX 



S \ TR 

Theorem II .11 is proved.! 



^y-^^y-^^daiy). 



\x - y\ 

s=\x-y\ 



3. Derivation of Theorems 11.21 and 11.3 



We prove Theorem 1 1.2 1 bv showing that it is equivalent to Theorem ll.il Indeed, 
it suffices to prove the following result: 

Lemma 3.1. Suppose that v G G[0,oo) such that v has proper decay at infinity, 
say v{X) = 0(|A|"^) as X oo. Then for any s > 0, 



(25) Gis,X)v{X)dX = Wiv)is). 



In order to prove this lemma, we need the following explicit formula for G: 
Proposition 3.2. Let c„ be as in For any s > 0, we have 



(26) G(s,A) = 



n-2 ( °^ 



\ s ds } 



, if n is odd. 

This formula must be well known. For completeness, its proof is given in Section 

m 

Proof of Lemma [O The relation ^ means that G(s,A) = W(e*^*)(s). 
Keeping this fact in mind, one sees that the following proof consists of just changing 
the order of the operator W with the integral sign: 
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• For even n, due to the decay of v and Proposition 13.21 the foUowing calcu- 
lations are valid: 

Gis,X)vWdX = c„/(i^) ' (J ^:L=dt \ viX)dX 

71 — 2 OO 



2 OO 

1 f v(t) 



• For an odd n, due to Proposition 13.21 and decay of v, 



G{s,X)v{X)dX = I [~] ( ) v{X)dX 



ld_ 

s ds 



1 d ^ f v{s) 
s ds 



v{X)dX 
Wiv){s). 



The lemma is proved. 1 



We now prove the following result, which shows that Theorem 11.11 implies The- 
orem [L3l 



Lemma 3.3. Let f G C^(IEl"), g = T{f) and go is its even extension with respect 
to t. Then 



00 

2 



G{s,X)go{yA)dX^- / A7^(s,A) / 7^s(/)(y, r)X(r, A)dr. 



We need the following auxiliary result: 
Proposition 3.4. Let f e C^(IR") and g = T{f). Then 

00 

g{y,s)e'^'ds = -tX J ns{f){y, s)G{s, X)dt. 


Proof: For the expository purpose, we provide here two proofs of this proposi- 
tion. 

00 

1) Let u(x, X) — J u{x, s)e^^''ds. Due to ([T]), u{., A) is the solution for 


AUix) + X^Uix) ^iXfix), 
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obtained by limiting absorption. Due to ([4]), we obtain 



u{yA) = -iA J f{x)G{\y - x\,X)dx 

oo 

= -iX f{nsf)iy,s)Gis,X)ds. 







This proves the proposition. 
2) Consider the transform 



i-^)^c,.J-s (i^)^ ( / ^dt) . if - i« even, 



( 1) ^ ^"ds il ds) ^ ( i '') 



if n is odd 



Let g £ C°°{Sx [0, oo)) , we also define B{g){y, s) = B{gy){s) where gy{s) - 
g{y,s). Let g — T{f), we then have (e.g., 013) 5 = BiJZsf)- Therefore, 



C50 OO 

g{y,s)e'^'ds = f B{nsf){y,s)e'^'ds. 



Substituting the expression of B into the integral and integrating by parts, 
we obtain 

oo oo 

iXs ^ 



g{y,s)e'^'ds = -iX J {Us f){y, s)W{ex)i.s)ds, 



where e\{s) — e^^" . Due to Proposition 13.21 we have 

oo oo 

J giy,s)e'^'ds = -iX J {nsf){y,s)G{s,X)ds. 



This completes the proof. I 

Proof of Lemma 13.31 Due to (|26l). G(s,— A) — G{s,X). Since 50 is even, we 
have 



(27) J G(s, A).9o(2/, X)dX = 2 J 7^(s, X)go{y, X)dX. 

u 

Since 170 is the even extension of g with respect to t, 

00 

5o(y,A) =2^/ 9o{y,t)e'^'dt^^Re J g{y,t)e'^'dt. 

JR 
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Due to Proposition 13. 4|, we get 



5o(y,A) = ^Re\^-iX I ns{f)iy,s)Gis,X)dt 

CO 

^ • ns{f){y,t)I{X,t)dt. 



TT 







From (P7)) . we arrive at 



G{s,X)go{y,X)dX^- f A7^(s,A) / 7^s(/)(y, r)X(r, A)dr. 



TT 

IR 



The proof is completed. I 



4. Proof of Theorem IL4J 



In this section, we present two proofs of Theorem II. 41 One of them rehes on some 
results of [TDl and works in the time domain while the other one is self-contained 
and works in the frequency domain. Without loss of generality, we can assume that 
n = B, the open unit ball. We need to prove that if / € C°°{B) and g ~ T{f) 
then for all x S -B 

(28) j jGi\x-ylX)go{y,X)dXda{y)^0, 

S H 

where go is the even extension of g with respect to t. 

4.1. The indirect proof. This proof is somewhat indirect since it uses inversion 
formulas in [lOl [9] . 

Due to Lemma [331 equality ([28]) is equivalent to 

W{go){yAx~y\)da{y)=0, 



s 



for aU X e B. Or, since g = 5o|sx[o.oo), 

(29) [w{g){y,\x-y\)da{y)=Q. 



s 

We introduce the following definition 

Definition 4.1. Let C{S x [0, oo)) be the space of all functions h e C°°{S x [0, oo)) 
satisfying the following conditions: 

i) h vanishes at t — to infinite order. 

ii) // n is odd then h is compactly supported. If n is even then for any non- 
negative integer k, \\d^h{.,t)\\[^^(^s'f — 0{t^"^''^^) as t ^ oo. 

Then the operator V, introduced in Section[l] maps C^{B) into C{S x [0,oo)). 
Indeed, one can show it by using the Kirchhoff-Poisson solution formulas for wave 
equation (e.g., [3). 

We now prove the following auxiliary result: 
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Lemma 4.2. Let x <E B and h e C{S x [0,oo)). Then 

V*{h){x)^ I W[h)[y,\x-v\)d<j{v), 



where V is the L^-adjoint ofV. 

Proof For an even n, a direct calculation (see [13]) gives 



oo 



1-2 



(30) V*{h){x) = c,, I I Ji^=^k^dtdy. 



y/t^ - \x-y\^ 



We observe that 



Id f h{y,t) ^ 1 d f h{y,t) t 



sdsj - s2 sdsj t y/P-s- 

S S 

OO 

= -l± fl(l^]VW^-dt 

sds J dt \ t J 

S 

OO 



d fh{y,t)\ 1 



dt\ t J 



Hence, by induction, for all fc > 0, 



HI) hiy,t) fldV f h{y,t) 



y/t^ - s2 \sdsj J yfW^ 



s 



2 



Therefore, due to (PU)) . we have 



-2 OO 



V*mx) . c„ / ( / -^$^dtdy 



sdsj J Vt^ - s2 ■ 



For an odd n, a direct calculation (see [TD]) gives 



s 

W{h){y,\x-y\)daiy). 



da{y) 

t=\y-x\ 



This proves the lemma.l 

This lemma tells us that is equivalent to V*{g) — 0. Equivalently, 
(31) V*dtV{f) = 0. 
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We now derive this equality from inversion formula in [TOl [9]. Indeed, if n is even, 
[21 Theorem 4] gives 

fix) = -2{V*td!Vf){x), 
fix) = -2iV*dttdtVf)ix). 

By subtracting these two inversion formulas we obtain the desired equality pip . If 
n is odd, [TOl Theorem 3] gives 

fix) = -2{V*dltVf)ix), 
fix) = -2iV*dttdtVf)ix). 

Again, subtracting these two inversion formulas, we obtain the desired equality 
(|3T|) . The proof of Theorem 11.41 is completed. 



Remark 4.3. As shown in [12 , identity I131\) is the complete range description for 
the operator V ( or, equivalently, for T ) when n is odd. 

4.2. The direct proof. Let ui be the extension of u by zero for i < 0, mi be its 
Fourier transform, and gi — mi|sxir,- Applying Lemma |3. 11 we have 

Gi\x ~ ylX)gQiy,\)d\d(jiy) = j Wigo)iy,\x - y\)daiy), 

S R S 

and 

j jGi\x-ylX)giiy,X)dXdaiy)^ J Wigi)iy,\x - y\)daiy). 

s n s 
Since go ^ gi on S x [0,oo), one has Wigo)iy, \x - y\) = Wigi)iy, \x - y\). The 
above two equalities then give 

i32)J jGi\x-y\,X)goiy,X)dXdaiy)^ J J Gi\x - y\, X)giiy, X)dXdaiy). 

5 R S R 

Due to Lemma 15751 

oo oc 

h{y, A) = ^ y giy, s)e'^'ds =^J 7^s(/)(2/, s)G(s, X)ds 



(33) I fi^)Gi\z-y\,X)dz. 

R" 

From ^ and we get 

Gi\x - y\,X)goiy, X)dXdaiy) 



s R 

^ J jGi\x^y\,X)X J Gi\y-z\,X)fiz)dzdXdaiy) 

S R R" 

^ I fiz) { I X I Gi\x-y\,X)Gi\y~z\,X)daiy)dX \ dz. 
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That is, 



(34) 1 1 G{\x ~ y\,X)go{y, X)dXda{y) }{z) j j X^{X)dX j dz, 

S M B VR / 

where 

(^(A) = / G{\x-y\,X)G{\y-z\,X)d<j{y). 



From the expUcit formula for G in Proposition l3.2[ we see that G(s, — A) — G{s, A). 
Thus, 

(35) ^(-A)=^(A). 

Let 7^(s,A) = Re{G{s,X)),J{s,X) = Im{G{s,X)), and 



if(a;, z, A) - J ni\x - y\,X)Ii\z - y\,X)daiy). 

s 

We claim the following symmetry whose proof is presented in Section [51 

(36) K{x, z, A) = K{z, X, A) for any x, z ^ B and A € H. 

Assuming this symmetry, one has Im{ip{X)) = K{x,z,X) — K{z,x,X) = 0. From 
we get f{—X) = (p{X). This implies the following equality: 

A^(A)dA = 0. 



Due to (IMl). one concludes that 



s u 



The proof is completed. 



Gi\x-ylX)go{y,X)dXdaiy)^0. 



5. Some special cases 



In this section, we derive Propositions 11.61 11.71 11.91 and ll.lOl from Corollary 1 1.51 
by proper choices of ^ and ip. 

5.1. Propositions [lT6l By choosing ^ — x and fix) = c in (fT2|) . we obtain 

fix) = -2 1 a.W(5)(y, s)L^|,_^| \x - y\da{y) + c J W{g){y, \x - y\)d<j{y). 

s s 
Here we use the notation dg for That is, 

(37) fix) = -2 / [sdsWg] (y, \x - y\)da{y) + c f W{g){y, \x - y\)daiy). 



We now claim an identity, whose proof can be found in Section [51 
(38) sdsWig) = Wisdsg) - (n - 2)^(5). 
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Assuming this claim, we see that ([BT)) is equivalent to 

fix) = -2 jw{sd,g){y,\x-y\)da{y) 

s 

+ [2{n -2)+c]J W{g){y, \x - y\)daiy). 
s 

Due to Lemma 1121 we obtain 

(39) fix) = -2V* (sd^g) (x) + [2{n - 2) + c]V*{g)ix). 

• If n = — 2(n — 2) then ([55)1 becomes 

fix) = -2V* isdsg) ix) = -2iV*sd^,Vf)ix) 

• If n = — 2(n — 1) then ((39| becomes 

fix) = ~2V* isdsg) ix) - 2V*ig)ix) = -2iV*dssdsVf)ix) 

• If n = — 2n then (|39l) becomes 

fix) = -2V* isd^g) ix) - iV*ig)ix) = -2(P*a^P/)(a;). 
Propositions 1 1 . 61 is proved (in this proof we have used the variable s in place of t). 
5.2. Proposition [T?fl Choosing fix) = in pT|) . we obtain 

(40) fix) f (^K^iy, s)) ^y-^^y-^^ daiy), 
Due to dSl) we have Gis, X) = ^"',1, [j J(As) - A^(As)]. Thus, 

4(27r)T~ 

(41) 7^(s, A) = - „_, A^(Ag), X(r, A) = —-—^ J(As). 

4(27r)— 4(27r) — 

Hence, 

oo oo 

K„iy,s) = J XnisX) jTZsif)iy,r)JirX)drdX 





oo 

-1 



16(27r; 



n-2 



A2"-3^(sA) J nsif)iy,r)JirX)drdX. 



Since / is supported inside B, we have TZsif)iy,r) = for r > 2. Therefore, 

cc. 2 

^"(2/'^) = / A2"-'A^(sA) / 7^s(/)(y,r)J(rA)drdA 



16(27r)"-2 
-1 



16(27r)"-2 
where, as defined in Section [U 

oo 2 

Kiy,s)^ f A2"-37V(sA) / 7^s(y, r)J(r A)drdA. 
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From (|40p. we arrive at 



s 



Choosing ^ — x, we have 



2(27r) 

That is, 



(43) /(X) ^ ^/(4.„„,,))_^^^_^^<i„(„). 

We claim that for even n 

(44) A;„(y,s)==i^i& /f-J^i) 7^5(/)(2/, r) In jr^ - s^jd^. 

TT J \dr r J 


Here the integral is understood in the principal value sense. A proof of this claim 
will be given in Section [6.41 Assuming it, we obtain 

, ^ 2(-l)^ f f d lY'\ ,,,, , , 



2 



2(-l)^ f f d V ""^ 



TT J \ dr r 





7^s(/)(y,r) 



dr. 



Since 

2 



we get 



7^5(/)(y,r)dr=- -- 7^5(/)(2/,r) 
dr r I r \dr ~ 



= 0, 





2 1 

r 





2 



2(-l)- /■ / d 1 



TT \dr r 







Integrating by part (which can be justified for the principal value integral in ques- 
tion), we obtain 



/ J 1 \ n— 1 



s#fc„(y,s) = 7^5(/)(y,r)ln|r2-.2|rf,. 



ds TT j dr \dr 



■ r 
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From (j43|) . we arrive at 

2 



/(^) 



(-1)- 



i2n) 



s 



dr \ dr 



d 1 



Equivalently, /(x) is equal to 

2 

d 1 



(-1) ^ UJ„ 
(27r)" 



s 



dr V rf?" 



n-l 



{y, r) In — — j/p drda{y). 



{y, r) In — \x — y\^ I drda{y). 



5.3. Propositions [T79l and [l.lOl Choosing ^ = in (|42p . we obtain 



/(^) 



2(27r) 



5=|k-;/| F 2/1 



-da{y), 



That is, 
(45) 



/(a;) = 2(27r)"-i ^"' 7 "(2/)fc«(y, |a; - 2/|)dCT(y). 



This proves Proposition ll.9l In order to prove Proposition ll.lOl we need only show 
that h(i/, s) = 2fc„(y, s) for any s > 0. Indeed, it is a consequence of the following 
lemma whose proof can be found in Section [6l 

Lemma 5.1. Suppose that h e C°°[0, oo) such that h does not grow too fast at 
infinity and /i^'^(O) = for all < i < n — 3 if n is odd, and /i^*^(0) = for all 
< i < n — 2 if n is even. Then 

OO C30 OO CXD 

X^"~^N{sX) J h{r)J{rX)drd\:^ ~ J X^'^-^JisX) J h{r)N{rX)dX. 



6. Proofs of auxiliary statements 
6.1. Proof of Proposition [3.21 We now prove 



(46) G(s,A) 



(7 T^dt], if n is 



Indeed, from ([3]), we get 



if n is odd. 



(47) 



Gis,X) -U:^) H^HiXs) = I \^_, . 
A\2TrsJ — 4(27r)— (As) — 



where H^^l2 is the Hankel function of the first kind. 

2 

• Let n be an even number. Using the equality (e.g. [25l page 74]) 



(48) 



1 dA™ /^ffi'^(.s)~ 
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we obtain 

H^'l (As) 



(-1)' 



(-l)^A-("-2) {-4- 
' s as 



Hence, (H71) gives 

4 (27r)^ Vsrfs 
Since (e.g., [H page 170]) 



is) ■ -»"<^) 



As 



i?(i)(As). 



U-1)^ 



1 d 



we conclude that 



-2i 



-2i 

TT 



- .S2 



G(s,A) 



(-1)' 



1 d 

s ds 



This confirms (j46p for even n. 

Let n be an odd number. Using again, one gets 



i/l'2.(A,s) 



= (-l)T^A-("^3) 



(As)' 

Hence, from (|T7|. we obtain 

G(s,A) = ^^^i^ 
4 (2^)^ 

Since (e.g., [3 page 487]) 



lrfX^ gf(A.) ' 
sc's/ (As)^ 

lrfX^ 4^(A.) 
s ds J (As)^ 



we conclude that 



G(s,A) 



1 d A ~ e*^" 
s ds / s 



2(27r) — 
This confirms for odd n. 
6.2. Proof of identity (|36p . We prove that for all A G IR, the function 



K(x, z. A) = y 7^(|x - y\,X)I{\z - y\,X)da{y), 
s 

is symmetric with respect to x, z € B. Here we follow the line of reasoning used in 

m- 

Due to (gll), we get 



(49) 



K{x,z,X)^ca / N{X\x-y\)J{X\z-y\)da{y), 
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where co is a constant depending only on n and A. We recall the following identities 
from [T8| for |a;| > tq: 



(50) 
(51) 



\v\=ro 



YiHy)JiX\x ~ y\)dy = Ci Jfe(Aro) Jfe(A|a;|)r,^(i), 
Yi'^{y)NiX\x~y\)dy - ci J,(Aro)7Vfc(A|x|)r,'=(x), 



\y\=ro 



where x — and YJ*^ is a spherical harmonic of order k. Consider the spherical 
harmonic expansion 

(k,l),(k' ,1') 

where a = \x\ and /3 = \z\. 
Due to (l49l). we obtain 



afc6'(a,/?) 



y,''(£)l^?'(i)i4:(af,/3z, A)dCT(i)dCT(z) 



s s 



= coy y Yi''{x)Yi';' (z) J N{X\y ~ ax\)J{X\y - f3z\)da{y)dxdz 
s s s 

^^ojij Yi\x)N{X\y-ax\)dA \ I YiY {z)J{X\y - (3z\)dz \ da{y). 
s \s 

Applying (fSOj) and ([5T|) . we arrive at 



a^;|,,(a,/3) = coc?Jfc(Aa)Jfe,(A/3)Jfc(A)7Vfc,(A)y r;:(y)rfc';(2/)da(y) 

s 

= coc2jfc(Aa)Jfc'(A/3)Jfe(A)iVfe/(A)4,fc'^M'- 

Hence, a^j(a,/3) = a^|;(/3, a) and K{x,z,X) = '}2(k,i)'h^\i'^^ PWt^^Wt^^)- These 
two equalities give K(x, z, X) — K{z, x, A). The proof is completed. 

6.3. Proof of identity (|38p . We now prove the identity 



(52) 

We first recall from ^ 



4w 

ds 



{g) 



w4 

ds 



{g) -{n- 2)W{g). 



!S^dt, if n is even, 



(53) yV[g){y,s) 



For n = 2, due to (l53l). we have 



if n is odd. 



oo 

4>V) (.)(., .)=c.(.i^)/ 



9iy,t) 
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Taking integration by parts, one gets 



giv^i) _ f d f g{y,t) 



- s2 J dt\ t 



Hence, 



oo 



''ds J ""J dt \ 



d ( g{y,t) 



-.dt 



dt\ t ) ^ 

oo 

^ ^Jl?^dt 



d fg{y,t)\ t^ r ^fgiy,t) 



dt \ t J - s2 y V ^ 

S S 

Taking integration by parts again, we arrive at 



d f g{y,t)__^^ ^ (dfgiy,t)\ e 



-At. 



ds J - s2 y V ^ J VP 

s s 

oo 



Simplifying the right hand side, we get 



d f 9iy,t)_j^^^ f tgt{y,t) 



dsj y/t^ - s2 J Vt2 

S 

Therefore, 



This confirms ([52]) for n — 2. 

Let n > 2 be even. We first observe the simple identity 
d \ ^1 d\ ./Id 



ds 1 \ s ds 1 \ s ds 1 \ ds J \ s ds 



By induction, we get 
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Therefore, due to the definition of W in (|53l) . 

n-2 OO 



S 



dds J \ ds J J \ft 



1-2 OO 



\ d\ ^ f g{y,t) 



-2 OO 



Due to (l54l). we obtain 



1-2 OO 



dds J J y/t 



= ( Ws- ) (5)(y,s)-(n-2)W(5)(2/,s). 



This confirms (|52p for any even n. 
Let n be odd. Due to (|53p. 



This and give 



d..\ . w . /i d\ W (g{y,s) 



ds J \s ds J \ ds J \ s 

71-3 



s ds / \ s 



Using the identity (sj^) (^'^^^^ = 5s (y, s) — ^Silifi for the first term of 
the right hand side of the above equahty, we obtain 



^-3 



ds 

This finishes the proof of identity ([52] 



= (ws4-){9){y,s)-{n-2)W{g){y,s). 
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6.4. Proof of (I44p . We prove a more general result: 
Lemma 6.1. Assume that n is even. Let g G C^[0, oo) and 

oo oo 

kn{y,s)= J X^"-^N{sX) J g{r)J{r\)drdX 



_T oo T 
(-1)^ 



Then 



kniy,s 



d 1 
dr r 



g{y, r) In \r — s \dr. 



Obviously, applying this lemma for g{s) — {TZs Div, s) we obtain (|44p. 
Proof: Since n is even, due to (pS)) . we have 



Thus, 



Since 



-— I iVo 



s ds 



(sA) / 5(j/,r) 



s ds 



--) Jo 



{r\)drdX. 



I d\~ 



and 



we have 



s ds 

I d\~ 

s ds 



Nn 



Jo 



(sA) = A-("-2) 
(sA) = A-("-2) 



s ds 



1 d\~ 



ds 



Jo(sA), 



kniy,s) 



1 d 



NoisX) / 5(2;, 



s ds 



1 d 
r dr 

1 d 



Jo{rX)drdX 



XNo{sX) I giy,r) ( -— ) MrX)drdX. 



Integrating by parts for the inner integral, we obtain 

(») 

where 



Uy,s) = / A7Vo(sA) / ( ) 5(y,r-)Jo(rA)drdA. 



Following the argument in [2], we obtain 



lniy,s) = -- 



d 1\ ~ 

dr ; 



9{y,r)- 



7 dr. 
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Here the integral is understood in the principal value sense. We now prove that 



+a 







for any nonnegative integer a. Indeed, we can rewrite 

OO 11-2 

2 f I / d 1 



r 



lu{y,s) = -[l — 9iy,r)— jdr. 

TT J r \ar r / _ 


Integrating by parts (which can be justified for the principal value integral in ques- 
tion), we get 

OO Ti-2 

lniy,s) = - / :^ — ( :^ — ) g{y,r)lii\r^ - s'^\dr. 
TT J ar r \ar r J 



Hence, 

OO 11 — 2 

s ds j ' TT I dr r \ dr r I ^ — 





OO 

1 [ ( d\ 

TT I \dr r 



9iy,r)^ ydr. 







By induction on a, we conclude that 



''n{y^s) = — / - — g{y,r)— — ^dr, 

s ds J TT J \dr r J r'^ — s'^ 



for any nonnegative integer a. Therefore, due to (j56p . 



Integrating by parts again, we obtain 

j^_2 OO n— 1 

k^{y^s)^tll^ J giy,r)ln\r^ - s^\dr. 





6.5. Proof of Lemma 15. 11 In view of (|4T|) . we have 

OO OO OO OO 

(57) J A2"-37v(sA) J h{r)J{r\)drd\ = c j >^nis,X) J h{r)I{r, X)drdX 



OO OO OO OO 

(58) J A2""V(sA) j h{r)N{rX)drdX = cj XI{s,X) J h{r)n{r, X)drdX. 
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Here the constant c is the same in these two formulas. Recah that, 

G{s,X)=W{e^)is), 
where ex{s) = e*^^. Let cosa(s) = cos(As) and sinA(5) = sin(As), we then have 

D OO OO OO 

A7^(s,A) J h{r)I{r, X)drd\ ^ J AW(cosa)(s) J h{r)W{smx){r)drdX. 



Similar to the argument in Lemma |3.1[ we can take W out of the integral sign to 
obtain 

OO C30 

(59) j A7^(s,A) J h{r)I{r,\)drd\ ^W{H){s), 



where 

OO OO 

H{s) = j Acos(As) j h{r)Wismx)drdX 



OO OO 

= J sin(As) J h{r)W{smx)drdX. 



Let W* be the L^-adjoint of W, which is 

, (-l)'^c„ f il^)^Lp^dt, if n is even, 

(60) W*{v){s)^l i ^ 

if n is odd. 



Given the condition in Lemma l5.ll integrating by parts (and changing the order of 
integration if n is even), we obtain 

OO 

h{r)W{sinx)drdX^ [ W* {h){r) smx{r)drdX. 



OO 

d 



Therefore, 



H{s) = sin(As) J W*{h){r)sm{Xr)drdX. 



Since the inversion of the Fourier-sine transform is itself, we arrive at 

ds 

Thus, due to 

OO OO 

(61) J A7^(s,A) J h{r)J{r,X)drdX = (^W^W*) Wis). 



Similarly, since 

OO OO OO OO 

/ A2:(s,A) / h{r)n{r, X)drdX ^ ( AW(sinA)(s) / h{r)W{cofix){r)drdX, 
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we obtain 



where 



AX(s,A) j h{r)'R.{r,\)drdX=W{Hi){s), 





Hi{s) 



Asin(As) J W*{h){r)cosiXr)drdX 



oo oo 

d 



ds 



cos(As) / W*{h){r)cos{Xr)drdX 



Here we have used the fact that the inversion of the Fourier-cosine transform is 
itself. Thus, 

CSO cso 

(62) j Xn.{s,X) J h{r)J{r,X)drdX = - (^W^W*) Wis). 



From dnZ]), dSSl), dH]), and we have 



X^n-iN(^sX) J h{r)J{rX)drdX ^ - J A^"-V(sA) J hir)N{rX)drdX. 



Lemma l5.ll is proved. 



7. Remarks 

The operator W is an intertwining operator between the second derivative 
and Bessel operator (e.g., [Ml EHl 13] ) , i-e. 



W 



n — 1 d 
s ds 



W(w). 



Thus, W transforms the wave equation into the Darboux equation, which 
is known to describe spherical means (see [6l [T4l l4l [15]V 
The symmetry of K{x, z. A) is similar to that of I{x, z, A) in [18| . It implies 
Theorem 11.41 as shown in Section 14.21 It is not clear that the converse 
is true. However, a close look at the argument in Section 14.21 shows that 
Theorem 11.41 implies the following weaker symmetry: 

J K{x,z,X)dX = J K(z,x,X)dX, for all a;, z G i?. 



R 



In jl3j , the authors derived an inversion formula using the known Neumann 
rather than Dirichlet observation data of the solution for wave solution for 
n — 3 (see [131 Theorem 12]). Applying the results in this paper, one 
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can derive the inversion formulas from Neumann data for any n. Indeed, 
looking at ([22]) and ((23|) . we arrive at: 

/(x) = 2 j jG{\y-x\,X)^iy,X)dXda{y). 

S TR ^ 

Due to Lemma |3. 11 we obtain 



fix) = 2 J W{d.^u)iyAx-y\)daiy), 
s 

where du stands for the (outward pointing) normal derivative. 
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